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Abstract 

With the first data sample collected by the CLEO-c detector at the ?/>(3770) resonance we have 
studied four exclusive semileptonic decays of the D° meson. Our results include the first observation 
and absolute branching fraction measurement for D° — > p~e + v e and improved measurements of 
the absolute branching fractions for D° decays to K~e + v e , ir~e + v e , and K*~e + v e . 
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The weak-current couplings of quarks within the Standard Model are described by the 
elements of the Cabibbo-Kobayashi-Maskawa (CKM) matrix 0], which must be determined 
experimentally. Because of their simplicity, semileptonic decays of hadrons provide powerful 
tools for probing the CKM matrix. Interpreting experimental measurements of semileptonic 
decay rates requires precision knowledge of form factors that are not easily calculated in 
quantum chromo dynamics because of non-perturbative effects. Form factor uncertainties 
are presently the main limitation in extracting \ V u b\ and \ V c b\ from semileptonic B decays 0. 

Heavy Quark Effective Theory (HQET) relates form factors in charm decays to those in 
bottom decays, and lattice gauge techniques calculate form factors in both charm and bot- 
tom decays. Precision measurements of semileptonic charm decay rates and form factors are 
a principal goal of the CLEO-c program at the Cornell Electron Storage Ring [3]. In this Let- 
ter, we report first results on semileptonic D° decays from CLEO-c: improved measurements 
of the branching fractions for D° — > K~e + v e , Ti^e + v e and K*~e + u e , and the first observation 
and branching fraction measurement for the decay D° — > p~e + v e . (Charge-conjugate modes 
are implied throughout this Letter.) The data sample used for these measurements consists 
of an integrated luminosity of 55.8 pb -1 at the ?/>(3770) resonance, and includes about 0.20 
million D°D° events jij. The same data and analysis technique are used for the branching 
fraction measurements of D + semileptonic decays in (Hf. 

The technique for this analysis, which was first applied by the Mark III collabora- 
tion at SPEAR, relies on the purity and kinematics of DD events produced at the 

'0(3770). We select events by reconstructing a D° meson in one of eight hadronic final 
states: K+tt', K+tt'tt , K+tt-jt ?! , K+ti-ji+ti-, K° s tc°, K^+tt', K° s tt + tT vr° , and K~K+. 
Within these tagged events, D° semileptonic decays are reconstructed in the exclusive fi- 
nal states: K~e + u e , 7r~e + z/ e , K*~e + u e , and p~e + u e , where K*~ — > K~tc° or Kgir~ , and 
p~ — ► 7r~7T°. Separation between signal and background from misidentified or missing parti- 
cles is achieved with the kinematic variable U = E m i BB — c\p m i BB \, where -E m i SS and p m iss are the 
missing energy and momentum of the D meson decaying semileptonically. The efficiency- 
corrected ratio of tagged events with semileptonic decays to the total number of tags gives 
the absolute branching fraction for the exclusive semileptonic decay mode. This branching 
fraction is independent of the luminosity of the data and benefits from the cancellation of 
many systematic uncertainties. 

The efficient reconstruction of tag events and the clean selection of semileptonic decays 
relies on the power of the CLEO-c detector, most components of which were developed for 
and used in B meson studies in the CLEO II and CLEO III experiments 0. The tracking 
system covers a solid angle of 93% of 47r with a six-layer low-mass stereo wire drift cham- 
ber surrounded by a 47-layer cylindrical (main) drift chamber. The main drift chamber 
provides specific-ionization idEjdx) measurements that discriminate between charged pi- 
ons and kaons. Additional hadron identification is provided by a Ring-Imaging Cherenkov 
(RICH) detector covering about 80% of 4ir. Identification of positrons and detection of neu- 
tral pions rely on an electromagnetic calorimeter consisting of 7800 cesium iodide crystals 
and covering 95% of An. 

Details of the criteria for selecting tracks, tt° and Kg candidates, and hadronic tags are 
provided in Ref. Q. The tag selection is based on two variables: AE = Ed — -Ebeam, the 
difference between the energy (Ed) of the fully reconstructed D° candidate and the beam 
energy (E^ eam ), and M bc = J 'E^ cam / c 4 — \pd\ 2 /c 2 , the beam-constrained mass of the D° 
candidate, where po is the measured momentum of the D° candidate. In case of multiple 
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TABLE I: The tag yields for the eight D decay modes with statistical uncertainties. 
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candidates, AE is used to select one D candidate for each tag mode, and we fit the beam- 
constrained mass distributions to obtain the tag yields. The signal component in these fits 
consists of a Gaussian and a bifurcated Gaussian to account for radiative and other effects. 
The background component is represented by an ARGUS function 8]. The yields of tags 
in all decay modes are given in Table |l] The total number of tags in our data sample is 
approximately 60,000. 

The requirement of a fully reconstructed D° meson tag greatly suppresses background. 
After a tag is identified, we search for a positron and a set of hadrons recoiling against the 
tag. (Only positrons are used because the CLEO-c muon identification system has poor 
acceptance in the momentum range characteristic of semileptonic D decays at the 0(3770).) 
Positron candidates are selected based on a likelihood ratio constructed from three inputs: 
the ratio of the energy deposited in the calorimeter to the measured momentum (E/p), 
dE/dx, and RICH information. Positron candidates are required to have momentum of 
at least 200 MeV/c and to satisfy the fiducial requirement |cos#| < 0.90, where 9 is the 
angle between the positron direction and the beam axis. The minimum momentum is 
chosen because of backgrounds from low-momentum pions. More than 80% of the positrons 
from D° semileptonic decays at the ■0(3770) resonance satisfy these requirements. The 
efficiency for positron identification has been measured primarily from radiative Bhabha 
events. For the criteria used in this analysis, it rises from ~ 50% at 200 MeV/c to 95% 
just above 300 MeV/c and then is roughly constant. The rates for misidentifying charged 
pions and kaons as positrons have been determined with exclusive hadronic decays of Kg 
and D mesons in CLEO-c data. Averaged over the full momentum range, the pion and kaon 
misidentification rates are approximately 0.1%. Bremsstrahlung photons are recovered by 
adding showers that are within 5° of the positron and are not matched to other particles. 

To select the hadronic daughters of a semileptonic D° decay, charged pions and kaons 
with momenta greater than 50 MeV/c are identified with criteria based on dE/dx and RICH 
information. Charged pion and kaon candidates must have dE/dx measurements within 
three standard deviations (3a) of the expected values. For tracks with momenta greater 
than 700 MeV/c, RICH information, if available, is combined with dE/dx. The efficiencies 
(~ 95% or higher) and misidentification rates (no more than a few per cent per track) are 
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determined with charged pion and kaon samples from fully reconstructed decays of D° and 
D + in CLEO-c data. 

We form tt° candidates with pairs of photons, each with an energy of at least 30 MeV 
and a shower shape consistent with that expected for a photon. The invariant mass of the 
photon pair must be within 3cr (er ~ 6 MeV/c 2 ) of the known 7r° mass. After selection, a 
mass constraint is imposed when ir° candidates are used in reconstructing other states. We 
reconstruct pairs of oppositely charged tracks from a common vertex to form a Kg candidate 



4.5<t) of the Kg mass. 



within 12 MeV/c 2 ( 

The K~ and 7r°, or Kg and tt~ , candidates are combined to form K*~ candidates. We 
require the invariant masses of the K*~ candidates to be within 100 MeV/c 2 of the mean 
K*~ mass. Likewise, 7r~ and 7i~° candidates are combined to form p~ candidates within 
150 MeV/c 2 of the mean p~ mass. 

A tag and the semileptonic decay are then combined. If there are no tracks other than 
the daughters of the tag and the semileptonic candidate in events, we compute U, which 
should peak at zero for a correctly reconstructed semileptonic decay. To improve the U 
resolution, we constrain the magnitude of the momentum of the D° candidate decaying 

semileptonically to be 



\/-^beam/ ° 2 



-f:0 



D° ' 



with its direction opposite to that of the tag D 
(p-po) in the i/>(3770) rest frame, where p-^o = — p D o. The distribution of U is approximately 
Gaussian with o ~ 10 MeV, varying mode by mode and somewhat larger for modes with 
neutral pions. In case of multiple K*~ or p~ candidates, we select a single combination based 



on resonance and 7r u masses. The U distributions for D° 
p~e + u e , summed over tag modes, are shown in Fig. ^ 
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FIG. 1: Fits to U = E m i ss — c|p m i ss | distributions for (a) D — > K~e + u e , (b) ir~e + v e: (c) 
K*~ (K~ir°)e + u e , (d) K*~ (K s ir~)e + v e , and (e) p~e + u e , with the D meson fully reconstructed. 
The solid line represents the total fit to the data, which includes the signal (dashed line), peaking 
background (dotted line), and non-peaking background (dot-dashed line). 



The signal yields are determined by fitting the U distributions. The signal is represented 
with a Gaussian and a Crystal Ball [9| function to accommodate the tails due to radiative 
effects. The tails of the signal function are fixed to the prediction of a GEANT-based ^(J 
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Monte Carlo (MC) simulation. The background function is determined from MC for each 
mode. The backgrounds are small and arise mostly from misreconstructed semileptonic 
decays with correctly reconstructed tags. The background shape parameters are fixed, while 
the background normalizations are allowed to float in all fits to the data. 

The background level for the decay D° — > K~e + v e (D° — > K*~e + v e ) is small. For 
D° — > K~e + v e , there is background from D° — > K*~ (fT~7r°)e + z/ e with an undetected tt°, 
but it is well separated from the signal because of the missing 7r°. For the decay D° — > ir~e + i/ e 
(D° — > p~e + z/ e ), the background peaks at positive U (well above the signal peak) and is 
primarily from D° — > K~e + v e (D° —>■ K*~ (K~7T°)e + i/ e ) where a kaon is misidentified as a 
pion. For D° — > K*~ (K~ 'ir )e + v e , there is considerable background from D° — > K^e + v e 
when a K~ is combined with a random tt° candidate, which is well below the signal peak. 
The signal yield for D° — > K~e + v e is determined by separately fitting subsamples for each 
tag mode. The results for all tag modes are found to be consistent. For the other modes, 
the yields are obtained with all tag modes combined due to limited statistics. The fits to the 
data are shown in Fig. and the yields are given in Table ITT1 The 31.1±6.3 D° — > p~e + v e 
events provide the first observation of this decay. 

The absolute branching fraction for any D° semileptonic decay mode is given by B = 
-^signal/ e-W"tagj where A^ signal is the number of D°D events with the tag D fully reconstructed 
and the D° reconstructed in that semileptonic mode, iV ta g is the number of D tags, and e 
is the effective efficiency for detecting the semileptonic decay in an event with an identified 
tag. Note that e = e si g na i/e tag is the ratio of the separate efficiencies for tag events with 
semileptonic decays and tag events in general. It is determined with a MC sample that 
includes relative populations of the eight tag modes that are consistent with the data. 
The cancellation of systematic uncertainties due to common effects in the numerator and 
denominator is clear. 

We consider the following sources of systematic uncertainty and give our estimates of their 
magnitudes in parentheses. The uncertainties in the efficiencies for finding tracks (0.7%) and 
for reconstructing it (2.0%) and K s (3.0%) are estimated with missing-mass techniques 0] 
applied to CLEO-c data and MC. The uncertainty in the positron-identification efficiency 
(1.0%) is taken from detailed comparisons of the detector response to positrons of radiative 
Bhabhas in data and MC. The effect of event complexity was incorporated by studying 
positrons both in isolation and embedded in hadronic events. The positron-identification 
efficiency depends on final-state radiation (FSR) and on bremsstrahlung in the material of 
the CLEO-c detector, the effects of both of which are simulated with MC. To assess the 
systematic uncertainty from these sources (0.6% combined), we vary the amount of FSR 
(simulated by PHOTOS [Tlj] ) and radiation in detector material, and we carry out the 
analyses with and without the recovery of radiated photons near positrons. Uncertainties in 
the charged pion and kaon identification efficiencies (0.3% per pion and 1.3% per kaon) are 
estimated by detailed comparisons of the detector response to tracks from hadronic D-meson 
decays in data and MC. There is an uncertainty in the number of 15° tags (0.7%), which is 
estimated by using alternative signal functions in the M^ c fits and by varying the end point 
(beam energy) of the ARGUS function parameterizing the background. The uncertainty 
in modeling the background shapes in the U fits (mode dependent: from 1.0% to 5.0%) 
has contributions from the simulation of the positron and hadron misidentification rates, 
as well as the input branching fractions. The uncertainty associated with the requirement 
of no extra tracks in a candidate event (0.5%) is estimated using fully reconstructed D°D° 
events in the data and MC. The uncertainty in the semileptonic reconstruction efficiencies 
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due to imperfect knowledge of the semileptonic form factors is small because of the uniform 
acceptance of the CLEO-c detector. It is estimated by varying the form factors in the MC 
within their uncertainties (1.0%) for all modes except D° — > p~e + v e , for which a conservative 
uncertainty (3.0%) is used in the absence of experimental information on the form factors 
in Cabibbo-suppressed pseudoscalar-to-vector transitions. The uncertainty associated with 
the simulation of initial state radiation (e + e~ — > DDj) is found to be negligible. Finally, 
there is systematic uncertainty due to limited MC statistics (0.7% to 1.5%, depending on 
mode). 

A non-resonant component is likely to contribute background in semileptonic decays to 
vector mesons. Based on evidence from the FOCUS experiment of a non-resonant component 
consistent with an S-wave interfering with D — > K*h>\ its contribution is estimated to 



be 2.4% in this analysis and subtracted when calculating the branching fractions for D° — > 
K*~e + v e . Systematic uncertainty associated with the subtraction (1.0%) is due to imperfect 
knowledge of the amplitude and phase of the non-resonant component. Interference with the 
S-wave amplitude alters the angular correlations among the decay products and introduces 
a systematic uncertainty (1.5%) in the reconstruction efficiency for D° — > K*~ (K~7[°)e + i> e . 
A relativisitic Breit-Wigner with a Blatt-Weisskopf form factor is used to simulate wide 
resonances in MC. A systematic uncertainty associated with the K* lineshape (1.2%) is 
estimated using D + -> K \K~ "7r + )e + z/ e , which has a much larger yield For D° — ► 
p~e + v e , there are insufficient data to constrain the non-resonant background or the resonance 
lineshape. Systematic uncertainties from these two sources are expected to be much smaller 
than the current statistical uncertainty for this mode, and they are neglected. 

The estimated systematic uncertainties are added in quadrature to obtain the total sys- 
tematic uncertainties in the branching fractions ( Table ITT]) : 2.8%, 2.9%, 4.9%, 5.1%, and 
6.6% for D° -> K-e+is e , D° -> n- e + u e , D° -> K*~ (K-n°)e + v e , D° -> K*~ (K%ir-)e + is e , 
and D° — > p~e + v e , respectively. Most of the estimates of systematic uncertainties are limited 
by data statistics and will be reduced with a larger data sample. 

In summary, we have presented absolute branching fraction measurements of D° semilep- 
tonic decays with the first 55.8 pb _1 of data collected with the CLEO-c detector at the 
^(3770). Our branching fractions for D° decays to K~e + i , e , 7r~e + z/ e , K*~e + i , e , and p~e + u e 
are given in Tables HQ and 11111 along with current world-average values compiled by the 
Particle Data Group (3. Recent CLEO III Q, BES and FOCUS |3 measurements, 
which are not included in Ref. ^3|, are also listed. The measurement of D° — > p~e + v e is 
the first observation of this mode. The absolute branching fraction measurements of other 
modes are more precise than and consistent with current world averages. Corresponding 
results for D + semileptonic decays and more extensive interpretation are presented in a 
companion Letter 0. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. This work was supported by the National Science Foun- 
dation and the U.S. Department of Energy. 
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